Supersymmetry representation of Bose-Einstein condensation of fermion pairs 
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We consider supersymmetry field theory with supercomponents being the square root of the Bose 
condensate density, the amplitude of its fluctuations and Grassmannian fields related to the Fermi 
particles density. The fermion number is demonstrated to be conserved in degenerated Fermi-Bose 
mixtures with unbroken supersymmetry when the system is invariant with respect to inversion 
of the time arrow. We show the supersymmetry breaking allows one to derive field equations 
describing behavior of real Bose-Fermi mixtures. Solution of related field equations reveals the 
cooling of homogeneously distributed fermions gives first spontaneous rise to strong inhomogeneous 
fiuctuations, while the Bose condensate appears at a lower temperature dependent of the fermion 
density. 

PACS numbers: ll.lO.-z, ll.30.Pb, 67.85.-d 



O 

o 



> 
o 
o 
cn 
p 

od 

o 

o\ 
o 



X 



I. INTRODUCTION 



Supersymmetry is one of the most beautiful and pro- 
ductive conception of contemporary physics which has 
been proposed for description of microworld (see Ref. [ij 
to review supersymmetry applications to the quantum 
mechanics, the survey deals with disorder metals, and 
book 3 relates to the theory of superstrings). Basing 
on idea proposed in the work [1] and developed in Refs. 
and @], it has been shown that the supersymmetry 
can be used effectively to present fiuctuative fields de- 
termining a picture of phase transition with symmetry 
breaking in non-equilibrium condensed matter 0] and 
random heteropolymers [1]. Along this line, the Martin- 
Siggia-Rose method ^ has been used as a basis per- 
mitting to combine stochastic fields into supersymmetry 
construction (lo| . One can perceive the supersymmetry 
approach is developed along two main lines: the former 
is based on the quantum operator representation [l|, Q , 
while the latter considers evolution of supersymmetric 

field 3i,i,S0,i,&li3- 

Initially, the supersymmetry had been proposed as a 
symmetry that relates bosons and fermions in elemen- 
tary particle physics [O, HI, [ll] ■ Along this hne, the su- 
persymmetric string theory is a unique theory expected 
to give a unified description of all interactions in nature 
Q , however none of the superpartners of any known ele- 
mentary particles has been found in experiments so far. 
Therefore, it is very important to study the supersym- 
metry breaking. Such an opportunity is opened with ex- 
perimental progress in atomic mixtures of ultracold Bose 
and Fermi atoms Theoretically, an ultracold super- 
string model was constructed jl5| . as well as an exactly 
soluable model of one-dimensional Bose-Fermi mixture 
was investigated [ll]. According to Ref. [l3|, the su- 



persymmetry is always broken either spontaneously or 
by a chemical potential difference between bosons and 
fermions. This article is devoted to consideration of the 
Bose-Fermi mixture within field-theoretical supersymme- 
try approach [T^ . 

The outline of the paper is as follows. In section [TTl 
we adduce main field-theoretical statements based on 
the generating functional method introduced by Martin, 
Siggia and Rose. Making use of this method allows us 
to write down both supersymmetric Lagrangian of the 
problem and Euler equations related. Section IIIII is de- 
voted to derivation of the field equations for the super- 
field components being the most probable values of the 
square root of the Bose condensate density, the ampli- 
tude of its fluctuations and the most probable Grassman- 
nian fields giving the Fermi particles density. Combining 
above equations, we show the fermion number is con- 
served in supersymmetrically degenerated systems whose 
superspace is invariant with respect to a rotation. Ac- 
cording to section IIVI the supersymmetry invariance is 
broken by means of fixation of such rotation that results 
in loss of the Grassmannian components invariance with 
respect to the time inversion due to the Bose-Einstein 
condensation. We show the supersymmetry breaking al- 
lows one to derive field equations describing behavior of 
real Bose-Fermi mixtures. Section|V]is devoted to discus- 
sion of obtained results and Appendix contains details of 
calculations at derivation of the field equations. 



II. MAIN FIELD-THEORETICAL 
STATEMENTS 

We consider attractive Fermi system characterised by 
the pair of conjugate wave functions ip{r,t), tp{r,t) and 
Bose condensate with the density n{r,t) where r, t are 
coordinate and time, respectively. Behavior of the Bose 
subsystem is presented by the fluctuating order parame- 
ter 
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with (j){r, t) being a condensate phase. Within framework 
of the standard field-theoretical scheme j^] , the system 
evolution is described by the Langevin equation 

dF 

x{r,t)^DV^x^~j — + ar,t). (2) 
ox 

Here, dot stands for the time derivative, V = d/dr, 
D and 7 are diffusion and kinetic coefficients, respec- 
tively, F{x) is specific free energy of condensate, C(r, t) 
is stochastic addition defined by the white noise condi- 
tions 

(C(r, t)} = 0, (C(r, t)C(0, 0)) = jTSir)5it) (3) 

where angle brackets notice averaging over system states 
scattered with variance T being the temperature mea- 
sured in energy units. With introduction of the scales 
= {-iTf/D\ rs = -/T/D, Fs = D^j^T^, (s = 
D^/{'yT)'^ for time t, coordinate r, specific free energy 
F and stochastic force C, respectively, the equation of 
motion ^ takes the simplest form 



xir,t) 



5T_ 

6x 



C(r,t) 



where short denotion of the variotinal derivative 
dT{x{r,t)} dF{x) 



6x Sx{r, t) dx 
is used for the Ginzburg-Landau model 

1 



T{x} 



F{x) + -(Vx)" 



V^x 



dr. 



(4) 



(5) 



(6) 



Along the standard line , our approach is stated on 
the generating functional 



Z{u[v,t)} 



/zWexp(/..drdt)D. 



being determined by the partition functional 



Z{x{v,t)} 




6F 
Sx 



CUet 



Sx 



(7) 



(8) 



where argument of the (5-functional takes into account the 
equation of motion ^ and determinant is the Jacobian 
of the transition from the noise field C(r,i) to the order 
parameter x{r, t) over whose distribution the continuous 
integration in the definition ([7]) is carried out. 

Within the simplest case of the Ito calculus, the Ja- 
cobian determinant equals one and the expression ([5]) 
arrives at the pair of the Bose field only Q . Much more 
interesting situation is generated by the Stratonovich cal- 
culus when the Jacobian 



det 



Sx 



exp 



d^^ = dV' (9) 



is presented by Grassmannian conjugate fields '0(r, t) and 
4>{j^, t) which subject to the following conditions: 



tptp + tp^p ^ 0, / dilj = 



dV- = 1, / dV' = 0, / V dV" = 1 



(10) 



Then, after generalized Laplace transform of the S- 
functional in Eq.® we obtain the supersymmetry La- 
grangian 



£{x,p, -0, '0) = ( pi - y + 

, , 9 S^F\ , 



(11) 



with a ghost field p(r , t) . Introducing the four-component 
superfield 



^ -.^ X + e^/j + iIjO + eop, 



(12) 



easily to convince that expression (|lip can be written in 
the canonical supersymmetric form 



A($) = ^(P$) (X>$) -1- F($) 



(13) 



where 0, 9 are Grassmannian conjugate coordinates de- 
fined by the properties 



ed0 = 1, de = o, ede = i 



(14) 



being similar to Eqs. ljlOp . Supersymmetry generators in 
Eq.([T3]) are as follows: 

d d d 

V:^Z-2e^, V:^-^. (15) 



III. SUPERSYMMETRY FIELD EQUATIONS 

Variation of the action related to the supersymmet- 
ric Lagrangian over the superfield (fT2|) derives the 
supersymmetry Euler equation 



1 r— 1 ST 



(16) 



where the square brackets notice commutation. As is 
shown in Appendix, projection of Eq. (jl6p onto the su- 
perspace axes 1, 9, 9 and 99 arrives at the explicit form 
of the equations of motion: 



r) — V f] 



8F 
8rj 



(17) 
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tp - 



drf 



Oaf 



(18) 



(19) 



(20) 



Since minimal action relates to the most probable realiza- 
tions of the superfield ([12]), solutions of equations p?)l - 
([20]) determine the most probable components = 

— max) = ^ j^j^^ :^('"'^^) ^ rpj^g 



7? 



first of these equations takes the form of the Langevin 
equation ^ that shows a ghost field ip = rep- 
resents the most probable realization of the fluctuation 
amplitude Specific peculiarity of the field cp is that 
gradient term in the governing equation (|18p has inverse 
sign, so that inhomogeneity in the space distribution of 
the most probable fluctuation increases in the course of 
the time until non-linearity stabilizes its amplitude. 

Another feature consists in the Grassmannian conju- 
gation of the equations (fTO)) and (pOI which coincide if 
the time arrow is inverted in one of them. Thus, one can 
conclude that the pair of conjugate fields ^ and ^E" de- 
scribes evolution of the Fermi particle and antiparticle for 
which the time runs in opposite directions. Combining 
Eqs. and ((20)) arrives at the continuity equation 



P + V j = 



for the fermion density 



and the current related 



(21) 



(22) 



(23) 



The equation (PT|) expresses the conservation law of the 
Fermi particles number in supersymmetric system whose 
state space spanned onto axes 1, 0, 9 and 96 is invariant 
with respect to a direction choice. 



IV. SUPERSYMMETRY BREAKING 

To break above invariance we shall follow to the Bo- 
golyubov method of quasi-averages, according to which 
taking off a system degeneration is provided by switch- 
ing an infinitesimal source type of slight magnetic field in 
magnets fl§\. In our case, the role of such a field is played 
by the conjugate Grassmannian components 9^ and ^'0 
the superfield (|12p which relate to forward and backward 
directions of the time arrow. Formally, we should replace 
the superfield (fT^ by the transformed field 



(24) 



Writing this superfield in the explicit form 

^ = X + {1 - x)^ + {1 + x)'^9 + {p + xp) 99, (25) 



easily to see that transformation ((24|) squeezes the axis 
6*5' and stretches the axis '^^9 by the same value x be- 
ing the order parameter of the Bose condensate, while 
the axis 99 is stretched by the value xp proportional to 
both order parameter and density of the Fermi particles 
p ~ ^"5. In any case, above transformation breaks the 
supersymmetry, so that the Euler superequation (|16p is 
reduced to the following components (see Appendix): 



dv ^2 dF 



dt 



(p + ryp) + VV 



d^F 
drf 



{p + tip) 



d^F 
drf 



(26) 



(27) 



vI'^ln[(l-77)vI/]_v2vI/ = -0vl/, (28) 



-^|ln[(l 



77) *] 



d'^F- 



drf- 



(29) 



In contrast to Eqs. (fT9|) and ([20]) . the pair of equations 
(PS)) and (P5|) becomes non-invariant with respect to the 
time inversion due to the Bosc-Einstein condensate ap- 
pearance (?7 ^ 0). Combining the equations ([28]) and 
(l29l). one obtains 



v") p] + Vj = 0. 



(30) 



At steady-state condensation (j 
at the relation 



: const). Eg. ([30)) arrives 



1-^ 



(31) 



where integration constant pc plays the role of a critical 
density of fermions and one takes into account the def- 
inition ([1]) according to which rf' = n. The dependence 
([3T|) means the density n of the Bose condensate increases 
steadily from n = to n = 1 with growth of the density 
p of Fermi particles above a critical value pc- 



V. DISCUSSION 

Characteristic feature of our consideration consists in 
making use of the supersymmetry field theory that is 
based on principle of the minimal superaction 

^{^(r,^} := j C[^{Y,t)]dvAt 
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whose values relate to maximal probability 

P{$(r, t)} oc exp (-5{$(r, t)}) 

in the system distribution over the superfields (fT^ . As a 
result, the governing equations ([^5]) - ([^5]) determine the 
most probable Bose components 77, (/? and Fermi ones 4", 
^. Such a description is differed crucially from the stan- 
dard picture, within whose framework observable values 
are determined in terms of averages over sets of quantum 
states. It is worthwhile to stress though the most proba- 
ble fields are determined without averaging over quantum 
fluctuations, however they take into account scattering 
over statistical states (see related averaging in Eqs.®). 

According to Eqs. and (PO)) . Fermi-Bose mixtures 
with unbroken supersymmetry are invariant with respect 
to inversion of the time arrow. To break this symmetry 
we transform the superfield to the form whose 
explicit appearance (j25p reveals breaking of above invari- 
ance due to the Bose-Einstein condensation. Because in 
macroscopical systems the time runs always forward, one 
needs to use the system ([^ - ([^ to describe real Fermi- 
Bose mixtures. 

To represent the system behavior let us study first the 
simplest case of steady-state homogeneous Fermi-Bose 
mixture. It is described by equations (pS)) . (|77|) and (I5D1) 
which are simplified to the forms: 



dF 
orj 



(32) 



The expression (j36p takes the familiar square-root form 

(38) 



Tc-T 
2T, 



cO 



in the Tc — T <C Tco vicinity of the critical temperature 

= (1 - ^/6^)j Teo. (39) 

Here, Eqs. ([5^ - give the stationary amplitude of 
fluctuations 

Ft tt" 

(40) 




and related density of fermions 



PO - Pc 1 + 



(41) 



According Eqs. p8|l and (|4T|) . with the temperature de- 
crease near the critical temperature (|39|) the densities of 
both fermions and bosons grow linearly, while the fluctu- 
ation amplitude (|40p takes negative magnitudes varying 
in square root manner. 

Let us consider finally inhomogeneous steady-state sys- 
tem that is described by the relation ((34)l together with 
the equations 



{e - ry2) 77 - (</? + r]p) 



(42) 



{if + 7?p) 



d^F 



V^) P, 



(1 -if) p = Pc 



(33) 



(34) 



To reveal analytically effect of external conditions we in- 
troduce the Landau free energy 



'2^ 



1 4 

4^ 



£ = 



Tco-T 



(35) 



with parameter e determining a moving off a character 
istic temperature Tcq. Then, solution of Eqs. ([32]) - ((34] 
arrives at the stationary order parameter 



2pc 



(36) 



It takes physically meaningful magnitude rjo ^ when 
the parameter e is more than the critical value 

ec = Voic (37) 

fixed by a critical density of fermions pc (together with 
the characteristic temperature Tco, above density pc rep- 
resents phenomenological parameter of the theory devel- 
oped). 



VV = - (e - 3r?2) {ip + 77P) - 6 (1 - 77^) r,p, (43) 



following from Eqs. and ^7} where the Landau 

free energy ([55)1 is used. Linearization of Eqs. (|1^ and 
(|43p over 77 and (p shows the e increase arrives initially 
(at e = 0) at loss of the homogeneity in space distribu- 
tion of the fluctuation amplitude and then (at e — Sc) 
the system becomes unstable with respect to the Bose- 
Einstein condensation. What about the fermion distri- 
bution, it is supposed to be homogeneous due to the 
equilibrium condition j = in the continuity equation 
([30)1 (consideration of more general steady-state condi- 
tion j = const 7^ is out of the scope of our study). From 
physical point of view, above means that with cooling of 
the Bose-Fermi mixture characterized by a homogeneous 
distribution of fermions strong inhomogeneous fluctua- 
tions appear spontaneously at a characteristic tempera- 
ture Tco, while following temperature decrease arrives at 
the Bose-Einstein condensation in the critical point Tc 
determined by Eq. (|39|) . 
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Appendix 



With accounting definitions (|T2)) and (|T5l). one obtains 
l.h.s. of Eq.lll 



(A.l) 



Being supersymmetry variational derivative, r.h.s. of this 
equation is written as the generahzation of the expression 



, , - — 
^ ' (A.2) 

- [v^Ty + e (v^*) + (v^f) e + ee (vV)] 

where the prime denotes the differentiation over related 
argument. According to the first rule (|14p. expansion in 
powers of the addition 0^ + "^9 + OQip gives 

F' ($) ^F'[ri+ (e^ +^e + OdLp)] 

= f' (r^) + f" (7^) {m +^e) (A.3) 
+ \F"{r])ip-F"'{r])p\ 69. 



Comparison of multipliers standing before 1,0,0 and 99 
arrives at the system of equations (fT7|l - ([20|l . 

In more complicated case of the transformed superfield 
(pS)) . the expressions (|A.1[) and (jA.SP take the forms: 



61* 



5- 

(1+7?) = + 77 



^'61 (A.4) 



at 



((^ + r/p) ( 



+ f" {f]) (1 - r/) + f" (r?) (1 + ??) *6' 

+ iv) + vp) - f'" (77) (1 - 772) p 



(A.5) 



Comparison of related supersymmetry terms derives to 
the system - 
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